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Abstract 

A theorem that constructs a path integral solution for general second order partial 
differential equations is specialized to obtain path integrals that are solutions of 
elliptic, parabolic, and hyperbolic linear second order partial differential equations 
with Dirichlet/Neumann boundary conditions. The construction is checked by eval- 
uating several known kernels for regions with planar and spherical boundaries. Some 
new calculational techniques are introduced. 
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1 Introduction 



The general path integral developed in [1] is specialized to yield path integral solutions 
to elliptic, parabolic, and hyperbolic linear second order PDEs. The examples in 
Section 3 provide explicit realizations of the construction and can be compared to 
known solutions. Some new calculational techniques are introduced which may offer 
some advantage in numerical methods. 



Material and notation presented in [1] will be assumed here. The reader requiring 
motivation to fully digest [1] may wish to begin with perhaps more familiar material 
in the examples presented here in Subsections 3.1 (the Poisson, diffusion/Schrodinger, 
and wave equations in unbounded space), 3.2.1 (the diffusion equation in the half- 
plane), and 3.3.1 (the Laplace/Poisson equations for a ball in R n ). As in [1] the 
issues of existence and uniqueness are not addressed; and functions, distributions, 
boundaries, etc. are generally assumed to be well-defined in any given case. 
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2 Path Integral Solution of PDEs 

2.1 General solution 



For reference purposes, the theorem of [1] is stated without proof. Relevant definitions 
and notation can be found in [1]. 

Theorem 2.1 Let Ml be a real(complex) m- dimensional (m > 2) paracompact dif- 
ferentiable manifold with a linear connection, and let U be a bounded orientable open 
region in M with boundary <9U. 1 Let f and cp be elements of the space of sections 
or section distributions of the (r, s)-tensor bundle over M. Assume given the func- 
tional S{x{r a i , z)) whose associated bilinear form Q satisfies Re(Q(x(r a r, z)) > for 

{x{T a ,,z))^0. 

//X(x a • E((r',r),z)) G F R (Q) where 

X(x a • r), z)) := / 9((r', r) - v)/(x a • E(r a >, z)) exp {-<S(x(r a /, z))} dr a , 

+ cp(x a • S((r', r), z)) exp {S(x((r', r), z))} ; (2.1) 
then, for x a = x(r a ) G U, 

I f x(x a -S(r^z))Vn (2.2) 
is a solution of the inhomogeneous PDE 

«F(x a ) = -/(x a ) (2.3) 

x=x a 

with boundary condition 

#(x B ) = <p(x B ) . (2.4) 

2.2 Elliptic PDEs 

Henceforth, restrict to the case where C + = M + or C + = iM.. Write r = st with s G 
{1, i} and f G IR+ or f G R which ever the case may be. For simplicity I will continue 
to write r in place of f with the understanding that r is real in this context. Although 
Theorem 2.1 holds more generally, it is appropriate to render the construction more 

1 The boundary <9U is assumed to be sufficiently regular. 
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accessible to standard applications and to connect with the notation of [6], [9], and 
[7]- 



Elliptic PDEs are characterized by closed boundaries with Dirichlet /Neumann bound- 
ary conditions and a positive definite Q. It is a simple matter to specialize Theorem 
2.1 for this case: 

Corollary 2.1 Given the hypotheses of Theorem 2.1, make two qualifications: Let U 
be a bounded orientable region in M with closed boundary <9U and assume the matrix 
G a(B has index (d,0). //x (£) (x a • E((t',t),z)) E F R {n) where 



X {E \x a -Z((T',T),z)) 



I 0(s[(t', r) - r a ,])/(x a • 27(7*, z)) expi^ 1 ^"^} dr a , 
Jc+ 



+cp(x a ■ r((r',r)^))exp{-^ 1 ^(^^^} ; 



then, for x a e U, 



!? (E) (x fl ) = / X {E) (Xa-Z(T^,z))Vn, 

is a solution of the inhomogeneous elliptic PDE 



-G^x (a) £x w + sC Y + V(x) 
with boundary condition 

^ E \x B ) = cf>(x B ) 



(2.5) 
(2.6) 



^)(x a ) = -/(x a ) (2.7) 



(2.8) 



The corollary clearly follows from Theorem 2.1. The only subtlety is keeping track 
of factors of s. The domain of integration C + is now R + or M depending on whether 
s — 1 or s — i. 

The kernels for Dirichlet/Neumann boundary conditions can be obtained directly 
from Section 3.2 of [1] and will not be repeated here. As an example, if the boundary 
of U is at infinity, V(x) — > V(x) + 2ir£ and = 0, then (2.6) can be written as the 
Fourier /Laplace transform of a path integral solution of an inhomogeneous parabolic 
PDE. The kernel of (2.7), for this case with vanishing Dirichlet boundary conditions, 
is 

K ( °\x a ,x a ,;£) = / / S(xa ■ £(T al ,z),x al )exp{-s- 1 S( y x(T a ,,z);£)}dT a ,V QjW z . 

JZ a JC+ 

(2-9) 

In quantum physics, (2.9) with s = i is the fixed-energy Green's function of the time- 
independent inhomogeneous Schrodinger equation when the boundary is at infinity. 
It is the Fourier/Laplace transform of the position-to-position transition amplitude 
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/C(x a /, r a r] x a , r a ) associated with the time-dependent Schrodinger equation. As previ- 
ously mentioned, the Fourier /Laplace transform interpretation becomes a Lagrange 
multiplier interpretation for phase space constructions. That is, (2.9) can be used ([7]) 
to give a phase space fixed-energy Green's function. 



Remark: In [7], the integrator Vr was chosen to be Gaussian, which in hindsight was 
not a good choice. However, comparing [7] with the present construction shows that 
the theorem in [7] still holds provided Dt is a gamma integrator. In light of Corollary 
2.1, the path integral in the theorem in [7] solves the inhomogeneous elliptic PDE 
with vanishing boundary conditions at infinity. Consequently, the fixed-energy Green's 
function calculated in [7] is the elementary kernel with Dirichlet boundary conditions 
at infinity and not the boundary kernel. 

2.3 Parabolic PDEs 

The path integral solution to parabolic PDEs when the (open) boundary <9U is 
nowhere tangent to the characteristic direction has been known for a long time. The 
solution, in the general context considered here, was presented in [6]. However, when 
a segment of <9U is tangent to the characteristic direction (for example, a boundary 
corresponding to a physical object whose position is fixed in time), supplemental 
Dirichlet /Neumann conditions are required along the segment. Both situations can 
be handled with some specialization of Theorem 2.1. 

2.3.1 General solution 

Because parabolic PDEs lend themselves to an evolutionary interpretation along the 
characteristic direction, it is convenient to assume that M is foliated so that the 
leaves coincide with a real (complex) (m — l)-dimensional paracompact differentiate 
manifold M~. The region of interest will be an orientable submanifold U which has 
the structure U~ xICM where U~ C M~, and the interval I C M + is a subset of 
the characteristic manifold of some parabolic PDE. The boundary of M is the union 
of two pieces; the boundary segment d\J~ xT := d\J~ x (I \ {0}) and the Cauchy 
surface U~ x {0}. 

It is necessary to construct a new parametrization for M. Let (x\ x°), % G {1, . . . , m — 
1} denote the coordinates in a local chart distinguished by the foliation. Instead of 
fixing the initial end-point, fix the final end-point; not because it is necessary but 
because it is instructive (and it conforms to physics usage). Hence, the final end-points 
of paths in M~ are fixed according to x{r b ) = x b e Mr and the final end-points of 
paths in I are fixed, x°(Tb) = x° 6 I. 
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Let T b + denote the space of L 2,1 functions r : T = [t a ,t&] — > [r a ,Tb] U [r a ,T5]* C C + 
such that \dr/dt\ > and r& := r(tb) = 0. Denote the space of L 2 ' 1 pointed paths z : 
r(T) -> M d (C d ) with z(t 6 ) = by Z 6 d . Let Q' = Zg x T+. Construct a parametrization 
P : Q' — > P 5 V M via the differential equations 



where the set of vector fields {X( a )(x(r(t))), Y(x(r(t)))} G T^^M - , the vector field 
Y°(x(r(t))) G T x ( r (t))I, and [X( a ) , Y°] = 0. Evidently there are two independent sets 
of parametrized paths; x{t, z) G M~ and x°(t) G I. Since the path x° will eventually 
be identified with the evolution parameter, take Y° = —d/dx® so the parametrization 
for x° according to (2.10) is just x°(t) = x° — t. 

Since the boundary has two pieces, there will be at least two critical paths. The 
critical paths x cr and x° cr reach dU at the first exit time t^. The critical path that 
reaches the Cauchy surface satisfies x° r (r^) = = x° — implying = x° for this 
critical path. Note that r^ b — > as x b — > x B or x° — > 0. 

It is now possible to construct the parabolic corollary to Theorem 2.1. 

Corollary 2.2 ^ssttme £/ie same hypotheses as in Theorem 2.1. Endow DCM with 
the foliated structure as described above and parametrize according to Eq.(2.10). The 
matrix G a/S is assumed to have index (d, 0). // x^Ki^b, x&) ■ E({t',t), z)) G ^(fi') 




o 

■6 



(2.10) 



( p )((x b ,x°)-r«r',T>,^)) 



:= / 9(s[(r',r) -r a ,])/((x 6 ,x°) •r(r a ^))exp{-^ l5 «^°^'^} aY a , 
+^((x fe ,x°) ■ r((r',r),,))exp{-- l5 «^°)(^^-))} , 



(2.11) 



£/ien ; /or (x 6 , x°) G U x I, 



(2.12) 



is a solution of the inhomogeneous parabolic PDE 
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.2 

Air 



S *G^C X{a) C Xw +sC Y + V(.)-s^ 



x=(x 6 ,x°) 



^ {P) (X6,X°) 

= -fM) 



with initial and boundary conditions 



x°=0 



v (x 6) o), v (p \x b ,x° b/ 



Xfc=X S 



<ax b ,x 



0\ 



(2.13) 



(2.14) 



Proof. It is not difficult to check that this corollary is a consequence of Theorem 
2.1. The only substantive difference is the presence of the d/dx b term in the partial 
differential operator, which is a consequence of the parametrization (2.10) and the 
choice Y° oc d/dx®. To show that the initial and boundary conditions follow from the 
theorem, it helps to use <9U = (<9U~ x IT) U (U~ x {0}) in condition (2.4), and to 
remember that — > as the boundary is approached. □ 



As a bonus, a "time-dependent" potential V(x a r,x®,) has been achieved. Also note 
that (2.13) is a diffusion-type equation for s real and a Schrodinger-type equation for 
s imaginary. 

2.3.2 Elementary kernels 

Lemma 2.1 The parabolic Dirichlet elementary kernel is given by 



K 



(D) 



((x 6 ,x°), (x a ,x°J:)^ J c+ 0(s[r^ - r a ,])<5(x 6 ■ Z{r a ,,z\ x a )5((x° - x°) - r a ,) 

xexp{-s~ 1 5((a;,a; )(v,2))} dr a ,Vtt' . (2.15) 



Proof. As in the elliptic case, the elementary kernel for the parabolic PDE is obtained 
by fixing the end-point with the delta functional 

6((x b ,x b )-S(r a/ ,z),(x a ,x° a )) . (2.16) 

From the parametrization for x° follows x°(r a ') = x° — r a i which implies the delta 
functional reduces to 

S(x b ■ Z(r a ,, z), x a )S((x° b - r„0 - x°) (2.17) 

with |x° — x°| > 0. It is straightforward to check that the integral has the correct 
form to apply the theorem, and that the kernel vanishes on the boundary. □ 
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It is useful to separate the boundary kernel into a kernel Kq V L on the boundary 
segment d\J~ x I~ and a kernel on the Cauchy surface U~ x {0}. 

Lemma 2.2 The parabolic Dirichlet boundary segment kernel is given by 

tf£?-((x 6 ,x°), (x B ,x° B )) := 8( Xb ■ E{r^ z),x B )5((4 - x° B ) - r^J 

xexpf-s-^^x )^^))} 2XT, (2.18) 
and the parabolic Dirichlet Cauchy kernel by 

^ D) ((x 6 , x°), (x a , 0)) := J S(x b ■ z), x a ) 

xexp{- S - 1 5((x,a; )(x°^))} £> Q ^ . (2.19) 

Proof. The boundary segment kernel clearly satisfies the homogeneous parabolic PDE 
and possesses the correct boundary condition. In specifying the Cauchy kernel, be 
mindful that the Cauchy kernel corresponds to a transition from a point x a on the 
Cauchy surface (x° = 0) to some point x& at the time x° — or vice versa: 

*f^((x b ,x°), (x a , 0)) := J w <5(x 6 ■ 27(7^, z),x a )5(4 ~ ^ ) 

xexp{- S - 1 5(( a :,x )(r 3 i^))} ZW' 

x exp {-^^((x, x )(r x ^, z))} d(\nr^ b ) V Q , w z 

= J zd S ( X b ■ ^( X 6>^)> X a) 

b 

xexp{- S -\S((:r,:r )(x°,z))} P Qi ^. (2.20) 

In the third equality, the constant M has been chosen to give K^q\{^, x°), (x a ,0)) 
the correct normalization. □ 

This is the well-known path integral kernel for parabolic PDEs for unbounded man- 
ifolds. For consistency, the Cauchy kernel must agree with the elementary kernel 
evaluated at x° = for the boundary at infinity. Putting x° = and — > oo in 
(2.15) indeed yields (2.19). 

Corollary 2.3 The solution to the inhomogeneous parabolic PDE with Dirichlet con- 
ditions 0(xs,x°) and Cauchy conditions i/>(x{,,0) is 
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# (P) (x 6 , 4) = I I 0(sK - (4 - x°)])<5(x 6 ■ E((4 - 4), z), x a )/(x a , x°) 

x exp x°)((x° - x°), ^))} VQ'dx a 

+ 8{x b ■ E{t^ z),x B ,)<f>(x B ,, (4 ~ r^)) 

xexp{-s~ 1 5((a;,a; )(r x ^,^))} VQ! dx B , 

+ / [ d 8(x b -S(x° b ,z),x a )^(x a ,0) 

J \J J z , 

b 

x exp ^-s' 1 S{(x,x Q )(xl, z))} V Q)W zdx a 

= j a 0(8[t£ - (xg - x°)])/(x b • 27((xg - x°), *),x°) 

x exp {-^((s, x°)((x° -*£),*))} W 
+ ^0(x 6 .r(ri^),(x°-r 5 i)) 

xexp{- S -\S((:r,:r )(Tx-^))} ^' 
+ J zd ^(x b -S(x° b ,z),0) 

b 

xexpj-s-^Or^Xx^))} V Q>w zdx a . (2.21) 

Proof. The expression follows from the two preceding lemmas. To check the boundary 
conditions, the first term vanishes for both x b — > x# and x° since — > in both 
cases. For the limit x& — > x#, the second term gives the boundary condition 4>(xb,4) 
along <9U~ x I - . The third term contributes i/)(xb,0) = 4>(xb,0) on dU~ x {0}. (In 
the third term, x b — > implies x° — > 0. This follows from the defining equation of 
Ifjf^x^xg), (xo, 0)) in (2.20).) When x° b -> 0, the second term gives no contribution 
because ^ 0. (It must be kept in mind that (x b — r^) > in the second term, 
because the integration over I~ did not include the point {0} G I.) The third term 
gives the initial condition ip(x a /,0). □ 

Observe that (2.21) satisfies the homogeneous PDE with Cauchy initial conditions 
only (with vanishing boundary conditions at infinity) by the particular choice f(x) = 
and <f){x) = 0. Specifically, 

«^ p )(x b ,x°)= J z y(x b -S(xlz))exp{-s- 1 S((x,x°)(4,z))} V QW z . 

b 

(2.22) 

With the identification x b = t b , this is the solution presented in [6] for the parabolic 
PDE with vanishing boundary conditions at infinity (as it should be). 
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The kernels for the case of Neumann boundary conditions for non-compact M can be 
adapted from Subsection 3.2.2 in [1]. 

Lemma 2.3 The parabolic Neumann elementary kernel is given by 

^r((x 6 ,x°), (x a ,x°J) xg), (x a ,x°)) + F u ((x 6 ,x°), (x a ,x°J) 



(2.23) 



where and Fu are now defined in an obvious way from (2.15). 
Lemma 2.4 The parabolic Neumann boundary segment kernel is given by 



' {N) ((x b y b ), (x B ,x° B )) := ^ 6(x(r^, z),x B )5((x° - x° B ) - r^) 

xexpl-s-'S^x^r^z))} TXt , (2.24) 

and the parabolic Neumann Cauchy kernel by 

tfg°((x 6 ,x°), ( Xa , 0)) := - 6(x b • 27(7^, *),x a )0(x° - r x \) 

xexp{- S - 1 5((a;,a; )(r x - L ft ,2))} Pfi' . (2.25) 

The proofs of these two lemmas will be omitted, because they are straightforward spe- 
cializations of the general Neumann kernels in section 3.2.2 of [1] and they essentially 
repeat previous arguments. 

Corollary 2.4 The solution to the inhomogeneous parabolic PDE with Neumann con- 
ditions V ns 0(x B ,x°) and Cauchy conditions V d / dx oi^(x b ,0) is 

# {P \mA)=I I ^{f ) ((x 6 ,x°),(x a ,x°))/(x a ,x°) 
Jv~ Jo.' 

x exp {-s^Sdx, x°)((4 - x°), z))} VQ! dx a 

+ L- L 6{x{t ^ z) ' xb ' )v "> (xb '' (x ° - 

x exp{-s" 1 5((x,a; )(r^,2))} VQ' dx B , 
~ j n , 0(4 - ri)V 9/dx o^ b ■ E{t^ z), 0) 

xexpj-s- 1 ^^ )^^))}^'. 

(2.26) 
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Proof. The terms follow from the preceding lemmas. It remains to check the boundary 
conditions for V na When x° — > 0, is in the x b direction. Neither nor 

Kq$- (since <jt 0) contribute. The third term gives the Cauchy initial condition 
VQ/Q x oi(y(xb, 0), because the derivative on the factor gives a delta function and the 

remaining 9 terms vanish for x b — > 0. When x& — > x^, V na «^ ( - p - ) gets no contribution 
from the first term and the second term gives the Neumann conditions V ng (f)(x B ,x b ) 
for x° > due to the derivative on the 6 factor and the vanishing of the remaining 
terms on the boundary. The third term contributes the Neumann conditions on 
the Cauchy surface V na V , ( x B ) 0) = V ns 0(x£,O) after recalling that x b — > xb implies 
x b — > (in the third term) and reducing the integral over T b + to a one-dimensional 
integral and integrating by parts. □ 

Note that it is possible to mix Dirichlet and Neumann conditions on the boundary 
segment and Cauchy surface. The expression for the associated «^ ( - p - ) will have the 
obvious collection of relevant kernels. 

2.4 Hyperbolic PDEs 

Like the elliptic and parabolic cases, hyperbolic PDEs can be solved by specializing 
Theorem 2.1. In fact, the solution is a rather trivial modification of Corollary 2.2. 

Assume that the matrix G a/3 has index (d—r, r) for some r e {1, . . . , d}. The quadratic 
form Q(x(r a r, z)) will have the same index since G a @ is ultimately determined by the 
inverse of Q. Recall that Re(Q) must be positive definite. It follows that the hyperbolic 
case will be associated with an imaginary Q. Furthermore, at least for wave equations, 
it is known that the boundary of interest must be open. (Whether this holds in the 
general case is, again, an open question.) With these qualifications, Theorem 2.1 can 
be applied directly. 

Because of its importance, it is useful to spell out the details for the case r = 1. It 
is convenient to formulate the corollary in terms of a foliated M as in Section 2.3 
and to fix the final end-points. The leaves of the foliation are assumed to be space- 
like. The evolution parameter is no longer along a characteristic direction, and the 
parametrization is determined by 

' dx°(r(t), z°) - Y°(x°(r(t), z°))dr = X°(x (r(t), z ))dz°(r(t)), x°(r b ) = x° 

< 

dx(r(t),z)-Y(x(r(t),z))dT = X (a) (x(T(t),z))dz a (r(t)), x(r b ) = x b 

(2.27) 

where the space of paths is Z£. 

Corollary 2.5 Assume the same hypotheses as Corollary 2.2 with the following qual- 
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ifications: The parametrization is determined by Eq. (2.27), the matrix G a/3 has index 
(d-1,1), and Re(Q(x(r a ,,z))) = 0. // x W ((*b, x°) • E{(? , t),z)) E F r (SY) with 
X iH H(* b ,4) ■ Z((t',t),z)) given by (2.11), then 



(2.28) 



is a solution of the inhomogeneous wave equation 



s 

An' 



: G^£ X(a) £x (/3) + sC Y + V(x) 
with initial conditions 



X=(X(,,X°) 



3^)(x b ,x°) = -/(x 6 ,x°) (2.29) 



Va/a x o^(x 6 ,0) , 



V d/dx o^ H \x b ,x° b ) 
and boundary conditions 



x o=o 



x fc =x B 



(2.30) 



(2.31) 



The kernels can be taken over from the parabolic case remembering that 

■ (w°J) * H** ■ £K,zU a )6((4 - 7^) - x°) 



(2.32) 



in general for the parametrization (2.27). So the Dirichlet kernels for the wave equation 
are just the parabolic Dirichlet kernels with 

5(x 6 -r(r x ^),x^(x°-^^ . (2.33) 

The Neumann Cauchy kernel is obtained by the substitution 

6(x b -IJ(T^z),xM4-ri)^S(x b -IJ(T^z),x^K),0) . (2.34) 

where 



0(x 2 ,xi) :-- 



(2.35) 



for x x 2< x 2 

1 for x x y x 2 
and the ordering is with respect to the foliation. 

The only task is to verify the correct boundary conditions. 

Corollary 2.6 The solution to the inhomogeneous wave equation with Cauchy data 
tf)(x b ,0) and J \/g/Q x oi/^(x b , 0) and Dirichlet boundary conditions <f)(x B ,x b ) is 
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* {H) M)= I /4 D) (K,^),(x fll x:))/(x ai x^x a ^ 

Jv Ji 

+\ [ if[f ) ((x 5 ,xO),(x a ,0))V(x a ,0)rfx a 
2 Jv- 

-i^K^ ) ((x 6 ,X°),(x a ,0))V W V(Xa,0) dx a . 

(2.36) 

Proof. To show the boundary conditions are satisfied, consider the limit x° — > in 
fl*^. The first and second terms vanish (recall that ^ I - ). The third term gives 
l/2-0(x6,O). In this limit, the fourth term can be integrated by parts (after reduc- 
ing the T fe + integral as before) to yield the remaining l/2'0(xfe, 0). The same limit in 
V d / dx o flK^ gets no contribution from the first two terms for the same reasons. The 
third term gives l/2V 9 / 9x oi/?(x 6 , 0) after reducing the T fe + integral and an integra- 
tion by parts. The fourth term contributes the remaining l/2V S / Sx o'j/>(x 6 , 0). When 
Xfc — > xg, the first term vanishes and the second term yields </>(xs,x°). The third and 
fourth terms give l/2if>(x B , 0) + l/2if>(x B , 0) = (p(x B ,0) as in previous cases. The 
fourth term requires an integration by parts as before. □ 

Neumann boundary conditions can be handled using the results of Subsection 2.3.2. 
The details will be omitted since they just repeat previous arguments. 



3 Examples 

As it stands, Theorem 2.1 should be characterized as a representation of a PDE 
solution rather than as a prescription for calculating the solution. The reason is that 
the complications introduced into the problem due to the geometry are encoded in 
r x , and one anticipates that techniques will have to be developed to handle these 
complications. Indeed, the parameter has been defined only implicitly. 

In the following examples, three such techniques are introduced in order to calculate 
the kernels for some selected geometries. The first is simply the observation that r x 
and T>t play no role when the boundary is at infinity. The second technique relies on 
the result of Appendix B to express the kernel of a bounded region in terms of the 
kernel on its covering space. The third technique is based on Subsection 3.2.1 in [1] 
and is akin to the method of images. 

In all of the examples, I will work in W 1 with Cartesian coordinates for simplicity. Since 
the target manifold is 1R™, it is natural to use the exponential map parametrization 
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which is explained in Appendix A. 2 Furthermore, I will take S(x) with V(x) = 0. 
Well known expansion techniques can be used for general V(x). 



3.1 Boundary at Infinity 



The easiest case to handle is when the boundary (excluding a possible Cauchy surface) 
is at infinity. Since the boundary is at infinity, — > oo (for any physically reasonable 
system at least). 



3.1.1 Elliptic Case 

For the elliptic case, take vanishing Dirichlet boundary conditions at infinity. Only the 
elementary kernel is needed in this case. Use the exponential map parametriza- 

tion for a classical path with Y = and X a = 5 l a -^ where % G {l,...,n} and 
a e {1, . . . , n}. 

The Dirichlet elementary kernel is given by 

K^\x a ,x a/ ) = J n J c S(x(r a ',z) -x a /)exp{-s _1 <S(z(T a >,2:))} dr a ,V^l . (3.1) 
The parametrization P yields 

x(r, z) = x a + z(t) . (3.2) 
The associated quadratic form is 

f T a' 

Q(x(r a >,z))= / 5 ij x l (T,z)x : >(T,z)dr 

J 

= r 5 afS z a (r)z^r)dr (3.3) 
Jo 

so that G a/3 = 5 a P. Consequently, K^ ) \x a ,x a i) satisfies 

^-AK { v D \x a x a ,) = -5(x a -x a ,) (3.4) 
where A is the Laplacian on R n . 

On the other hand, the exponential map parametrization gives 

x(tX) = x cr (r) + C(t) (3.5) 



2 Of course, the final expressions do not depend on the particular parametrization that one 
chooses. 
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where critical point of S relative to a variation with both end-points fixed, 



viz. 



x cr (r) = x a + ( — —Jr. (3.6) 

(Since the boundary is at infinity, this represents a transition between interior points 
so the relevant variational problem is for point-to-point.) 

The Gaussian integrator is invariant under translations, and so the integral can be 
expressed in terms of C, instead: 

K { ® } (x a ,x a /) = / / 8{x{r a i, C) - x a /) exp {-s -1 7rQ} dr a , Vtt , (3.7) 



where ' 



Q(x(r a 'X)) = Xa ' 2 + S af3 C a (r)C\r)dT . (3.8) 

To' Jo 



To evaluate use' 



S(x(T a ,, C) ~X af ) = 6(x cr (Ta>) + C(Ta') ~ *«>) = I 1 exp^ 2 ^^'^ du 

Jm. n 

= j Rn lexp{-2m(u5 C(Ta/) X)}du (3.9) 

to get 



K\?\x a ,x a ,) = J„J^ ^lexp{-s 1 ttQ - 27ri((u<J^ (v) ,C)} dudT a ,Vtt 

„ „ f ~ 7r l x a /~ x al 2 1 

=J T J c+ le{ STa ' ] L^{- s * w ^j) d ^ v - 

= [ [ lei '"v ' } / expH^M dudT a ,VT 

JT+ JC+ JR n 

r r f -^|x a ,-x a | 2 \ 

= / / lei "V J( S r a ,)- n/2 rfr a ^r 

„ f - 7r l x a '-*a| 2 1 

= / lei J( S r a 0- n/2 rfr a , . (3.10) 

For s = 1, this evaluates to 



3 There is no cross term since £(0) = for the exponential map parametrization, and 
C( T a') = will not contribute by virtue of (3.9). 

4 Recall from [1] that S(x a ,x a i) denotes a Dirac bitensor composed of Kronecker delta 
symbols (which are collectively denoted by 1 for brevity) and the scalar Dirac delta function 
<5(x a ,x a /). 
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-2 In \x a i — xJ 1 for n = 2 



(3.11) 



and for s — i, 

K { ®\x a ,x a ,) = 7(x a ,x a /;n)l (3.12) 

where 

/(x a ,x a ,;n) := r n/2 2 n/2_1 / k (n/2 ~ 2) exp |27ri k|x a , - x a | 2 } rfk (3.13) 

after the change of integration variable (r a /) _1 — > 2k with k real. In particular, for 

n — 4, 

J(x a ,x a ,;4) = 25[|x a , -x a | 2 ]l. (3.14) 

Note that the same answer obtains if the parametrization x(t, z) = x a + z(t) is used 
instead of the exponential map parametrization. In this case, the action functional 
has no |x a — x a /| dependence, but the delta functional does. Writing the delta func- 
tional as an integral and evaluating the resulting functional integral yields the correct 
exponential dependence on |x a — x a /|. 



3.1.2 Parabolic Case 

Here the boundary segment at infinity is <9U~ — not the cauchy surface. In this case, 
vanishing Dirichlet boundary conditions necessarily implies lim Xfc _ >00 i/^x;,, 0) = since 
-0(xb, 0) = 4>{xb, 0) must hold for consistency. The parametrization is the same as in 
the elliptic case with the addition of x°(t) = x° — t. The x° sector contributes nothing 
to Q since / \x° — Y°| 2 dr = 0. (This follows from the parametrization (2.10).) 

The kernel associated with transitions from the Cauchy surface is given by (2.19). In 
this example, it satisfies 



s 2 A d 

—A - s- 



4tt <9x° 
The calculation proceeds as before; 



K ( c D \(x b ,4),(x a ,0))=0. (3.15) 

X=(x a ,X S ) 



Kg ?) ((x 6 ,x°),(x a ,0))= / d d(x b -E(4,C)-Xa)e W {-s-\Q} V^C 

JT Xcr Z, 
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= K)-/ 2 exp ^V 6 ' [1 (3-16) 



sxi 



for x? > 0. 
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For the elementary kernel find, using (2.15), 



K{f ) ((x 6 ,x°),(x a ,x°))= If f 5(x b .i7(v,C)-x a )5(x°-x°-r a 

JT a JT Xcr Z b JC+ 

x exp {-s _1 7r(5} dT^V^^Dr 
= / / 16(4 - x° - T al ){sT al )- n ' 2 e\ /rfr^r 

"'To "^ C+ 



= / l5(4-x° a -T a ,)(sT a ,)- n/2 e\ Ur a , 

J C-L 



(3.17) 



For 8 = 1, l4 D) ((x 6 ,x°), (x a ,x°)) = for (x° - x°) < 0. If (x° - x°) > 0, then 



K^,4),(^X)) = (*£-x°r /2 exp{ ^jf Y ■ (3-18) 
For s = i, 

^ D) ((x 6 ,x°),(x a ,x°)) = [»(xg-xg)r/ 2 exp| ^J^ |2 |l (3.19) 

for either (xg — x°) > or (x° — x°) > 0. Intuitively, for imaginary s, — > -K"^ 
under time reversal. 



3.1.3 Hyperbolic Case 

Consider the wave equation. For this case Q is imaginary with index (n — 1, 1). Choose 
Y = Y° = 0, X° = <9/<9x°, X a = 5^-7^1 and the coordinates in the obvious way so 
that the PDE for the elementary kernel is 

*Ti D) ((x 6 ,x°), (x a ,x°)) = -<5((x 6 ,x°) - (x a ,x°)) (3.20) 

x=(x 6 ,x°) 

where 



47T 



A - 



d 2 



Kg\(*,,4),M))=LL s((* b ,4)- ^,0-^,4)) 

Jil J C7_|_ 

x exp dr a iVVl . 



(3.21) 



The quadratic form is dictated by the form of the PDE to be solved. It is more 
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compact to express it in terms of coordinates y p = (y°, y 1 , . . . , y n x ) on E( 1,n ^ with 
metric r\ p(J : 



~ f T a' 

Q(y(r a >, C); s) = j o VpaV P (r, z)y a (r, z)dr 



= + r V^ a {r)c\r)dr . (3.22) 

The first term in the second equality comes from the critical path in the same way 
as the elliptic case. 

Repeating the calculation for the elliptic case yields 



K ( v D \y b ,y a )= Iff lJy (y ^' b) hxpWv} dudr^Vr 

JT+ JC+ JK. n 

= 1 lei i(sr a/ )- n / 2 dr a/ Vr 
Jt+ Jc+ 

= / lei «v J( S r a 0- n/2 rfr a , 

/■ f -^'(ya-yi,) 2 \ 

= / lei v J(ir a 0- n/2 dr a , 

= /(y 6 ,y a ;n)l (3.23) 

where the second-to-last equality follows because Q = —Q* allows s — % only. In 
particular, in E^ 1 ' 3 ), 



K < v D \j b ,j a ) = 25[(j a - Jb ) 2 ]l. (3.24) 

For y° ' > y , 5 this corresponds to the retarded Green's function which has support 
on the future light cone. Since the kernel is symmetric under j a <-> y a >, the same 
kernel obtains for y° > y° ' and corresponds to the advanced Green's function. 



3.2 Planar Boundaries 



For planar boundaries, the issue of dealing with can be side-stepped by employing 
a trick. Express U as quotient space U = U/G where U is isomorphic to R™ and G 
is the discrete group required to obtain the desired boundary. It turns out that the 



5 Recall that either x° > x° or x° > x°. 
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elementary kernels on U can be expressed in terms of (presumably) known quantities 
on U. 

Let U denote the universal covering of U. The manifold U is a principal fiber bundle 
with projection II : U — > U = U/G and discrete structure group G. The parametriza- 
tion is constructed via the composition P = P u o P where P : Z% — > PjlJ and 

Given a connection on II and a fixed initial point x a G U, a given path x(-,z) can 
be horizontally lifted into the set of paths {x(p)(-, z)} with the index p running over 
the order of the discrete group. Hence, defining the evaluation map e : Z% — > U and 
its lift e : Z% — > U, then it follows that £ = II o e and e r_1 (x( p )(-, z)) G Z^ p where 
Z^p parametrizes the pth path x^(-,z). Consequently, the inverse image of the set 
of paths with a given base space end-point consists of the disjoint union of inverse 
images of the set of paths with end-points in the fiber over the base space end-point; 
i.e., Z^: can be decomposed as {J p Z^ p with Z% p fl = \/p ^ q. 

In Appendix B, it is shown that the elementary kernel on U can be expressed in terms 
of the elementary kernel on U (provided it is an equivariant map) by 

K r v (x a ,x a ,) =J2p r (9p 1 )K v{p) (x a ,x a ,) (3.25) 

v 

where the p r furnish unitary representations — labelled by r — of the structure group 
G, and -Ku(p) is the elementary kernel on II associated with the integral over Z% p and 
is given by 

Ku {p) (x a ,x al ) = / 9(r^ a - T a ,)5(x a ■ [g^E^a^z^gp]^) 

J £lp J C_|_ 

x exp {-S(x(T a > , z p ))} dr a , VQ P 

(3.26) 

where z p G Z% p and £l p := Z% p x T+. Apparently, using this trick and variations of the 
following three examples, one can find (at least in principle) the elementary kernels 
of fairly general n-polygons (notwithstanding smoothness issues on the boundary). 



3.2.1 The Half-Space 

In this example, I treat only the parabolic case. The elliptic and hyperbolic cases 
can be obtained similarly. Define the half-space U~ = W] := {x = (x 1 , ...,x n ) G 
M. n \ x n > 0} so that U~ = M. n . The boundary segment <9U~ is at infinity for points 
(xV..,x»-\0). 
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In this case, G = Z 2 and the parametrization P is given by 

*(t, C) = (xl r (r) + C\r), . . . , \x n cr (r) + C(r)|) . (3.27) 

Since G = Z 2 , the inverse image of a point x is x± := (x 1 ,...,±x n ) where the 
parametrization P gives 

±x n (;z) = ±(x%+z n (-)). (3.28) 
Defining x± := {x 1 , . . . , ±x n }, the delta functional contributes 

S(x a/± ±z(T a ) -Xa) . (3.29) 

Using the results of Subsection 3.1.2 and (3.25), the elementary kernel is 

^((x a ,,4),(x ,x°)) 

= P r (#+)^U(+)((Xa',X°,), (Xa,X°)) - p ? ' (#_ ) K V{ _) ( (x a > , X°, ) , (x a ,X°)) 

(3.30) 

where 

K u(±) ((x a ,,4),(x a ,x°J) = K4-^r" /2 exp( ~"j X r^f }l • 

(3.31) 

for real and imaginary s (with the appropriate restrictions on (x°, — x°)). 

For the trivial representation, [p^ D \g + ), p^ D \g_)) := (1,1), ifu vanishes on the 
boundary x™ = 0. For the representation, (p( N \g + ), p( N \g_)) := (1,-1), V n9 -K\j 
vanishes on the boundary x™ = 0, and, therefore, these two representations yield the 
Dirichlet and Neumann elementary kernels respectively. 

Alternatively the Dirichlet and Neumann kernels for the half-space can be obtained 
directly from Proposition 3.1 in [1]. The transformation a is just a translation of 
x™, — * so that x™ — > x™ + x™,. Making use of Subsection 3.1, one easily arrives at 
(3.30). 

The Cauchy kernel follows similarly; 



^((x 6 ,xj),(x a ,0)) 

= /(( ?+ )K c(+) ((x fe ,x°),(x a ,0))-/(^)K C ;H((x 6 ,x°),(x a ,0)) (3.32) 
where 
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K c(±) ((x 6 ,x"),(x a ,0)) = W)-^exp 



°W 2 — / z!!fe ~ Xfe± H i 



(3.33) 



The Dirichlet and Neumann kernels obtain for the appropriate representation of Z 2 . 
To find the boundary kernel 6 , calculate 



K 



(D) 

aw- 



((x b ,x°), (x B ,x° B )) = [6{(4,y$) ■ S{t^ z),x B )<5(x B - x° + r^) 

xexp{-s- 1 5((x,x°)(r 3 i^))} Ptf (3.34) 



where x B = (x B , 0) with i G {l,...,n — 1}. The boundary kernel vanishes for x B ^ 0. 
Since x B = 0, the integral over Z b reduces to the unbounded case. Moreover, the T b + 
integral can be reduced according to (B.54) in [1] since it depends only on r^; 



l^u4(x 6 ,x°),(x B ,x B )) = / r+ l«) 



±\-n/2 t 



-Trjxg-xi] 

IT 1 



x5(x° B - 

:AT( S (xO-X B ))'( n / 2+1 ) f 



-7r|x B -x 6 | 2 


} 


"A 








J -t|x S -x 6 




1 s(*r^> 





(3.35) 



for x° > x B and IfgJL 



otherwise. 

'(D) 



To determine A/", normalize K 9V L by choosing <fi(x B ,x° B ) = a where a is a constant. 
Then 



= aU I 6 l^xg-x"))- 3 /^!^™ / dx° B 
Jo 



erf < 












i 









= aAT(sx™)- 1 

Since «F (p) (x B , x°) = a and erf(0) = 0, then M = -sx™. Finally 



(3.36) 



6 Admittedly, the boundary kernel can be obtained from the normal derivative of 
but it is interesting to derive it from a path integral. 
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(3.37) 

for x° > x° B and Kg$_ = otherwise. 

More generally, this procedure can be applied to M n "folded" into l/(2m)-space (m < 

n) by the structure group Z 2 ©,...,©Z 2 where each "fold" is along a different x\ 

. ' 

m 

Another possibility is to use Z h with h a positive integer with U = R 2 to obtain a 
sector in the plane. 



3.2.2 The Unit Strip 

Consider only the elliptic elementary kernel case for brevity. Define the unit strip by 
U = {x = (x 1 , . . . ,x n ) e M n | < x™ < 1. An appropriate parametrization given the 
initial point x a = (x a , . . . , | sin(kx^) ]) =: (x a , . . . , x") is 

*(t, = (xM + C(r), . . . , | sin [k(x»(r) + C"(r))]|) (3.38) 

where x cr (r) is given by (3.6) and k is a real constant. The structure group is G = 
Z 2 © Z since the set of inverse images is {x(±,m)} = ( x \ • • • , ± x m) where 

±5* z) = ±(5C + ± 2vrm/k) (3.39) 

with m/k G Z. The delta functional is 

*(*o(±,m) ± ^( T a0 - x a') (3.40) 

with x (±>m ) := (x 1 , . . . , ±(x n ± 2-Km/k)). 

Repeating the relevant steps from Subsection 3.1.1, find the elementary kernel; 



K\J(±,m) ( x a,X a ') 

7rW2 r (2 _ i)| Xo , _x a(±)m) | 2 -n forn + 2 
-2 In |x a / - x a(±>m) 1 1 for n = 2 

= / ( x a(±, m ),x a ';n)l (3.41) 

for s = 1 and s = i respectively. Finally, for the trivial representation of Z, the 
Dirichlet and Neumann elementary kernels are 
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^Hj( x fl, x a') = [-K'u(+,m)( x a,X a ') =F K u(_ m ) (x a , X„/) 



(3.42) 



5.^.5 JTie E/ra* 5ox 

Let U be the box defined by U = {x = (x 1 , . . . , x n ) e R n | < x< < 1 Vi G {1, . . . ,ra}}. 
I will calculate only the elementary kernel for the elliptic case again for brevity. For 
the initial point x a := (| sin(kix*)|, . . . , | sin(k n x")|, the parametrization is 



x(r, C) = (| sin MxKt) + C l (r))}l . . . , | sin [k n (x n cr (r) + C(r 
with x cr {r) given by (3.6). The structure group is 

g = z 2 © z, e, . . . , ©z 2 © z , 



(3.43) 



(3.44) 



and the set of inverse images is 

{x(±,m)J = (ix^,...,©^) 

where 

±ai H (-,z)=±& + z i (-)±2irm i /*i) 
with i G {!,...,»} and rrij/kj e Z. 



(3.45) 
(3.46) 



Proceeding as in the previous examples, find the kernel for s = 1 and s = i respec- 
tively: 



= < 



n i-n/2 r( a _ 1))Xa , _ x a(±;m) .| 2 -1 for n ^ 2 



-21n|x a / -Xo( ±>m )Jl 



for n = 2 



: -T(xo(± ) m) j » x o';ra)l 



(3.47) 



where X( ±im ) j := {±(x 1 ±27rmi/ki), . . . , ±(x n ±27rm ri /k n )}. This yields the Dirichlet 
and Neumann elementary kernels (for trivial representations of the Z) 



K^(x a ,X a ,)= [■K'u(+ 1 m) j (Xo,X /) © ^(-^.(SCo.Xo/ 



(3.48) 



where {to 3 -} := {mi, . . . , m n }. 
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3. 3 Spherical Boundaries in W 



3.3.1 Then-Ball 

To handle spherical boundaries, it is necessary to come to grips with ■ The cal- 
culations for the boundary kernels can be simplified substantially by expanding the 
paths about a critical path and then using (3.33) in [1] which allows the point to 
be replaced by an arbitrary xb in the boundary kernel. 

Consider Laplace's equation for the case of the n-ball in M n defined by B™ := {x = 
(x 1 , . . . , x n ) G R n \x 2 — R 2 < 0}. First calculate lfg^(x a ,xf ). As usual, it is expedient 
to expand about a critical path and the exponential map parametrization is 

x 2 (rX) = [x cr (r) + C(r)] 2 (3.49) 

where 




(3.50) 



is a critical path. 

The form of x cr is dictated by the initial condition of the parametrization, which 
applies to x(t, £); and by the transversality condition (see Section 2 in [1]), which 
requires that cc cr (r^) intersects the boundary transversally. A reasonable choice is 
\&cr{ T x a ) I \ x a\/ T x a - Note, however, that the transversality condition does not fix the 
proportionality constant. The proportionality constant is ultimately fixed by requiring 
thatKf ) (x^,xf) = < 5(x B ,-xf). 

Now evaluate 

Kf\^ a ) = I \8{x{rtX) -x^expj-rtg} VQ , (3.51) 
where the quadratic form is 



Q(C ri) = / Sijx^r, z)x j (r, z)dr 

J 

= |X " ~ Xa ' 2 + 5 a ,C{r)c\r)dr . (3.52) 

7VT JO 



The first step is to note that the integrand depends on only. Therefore, it reduces 
to an integral over C + as usual. Using (B.54) in [1], get 



Kf\^ a )=N f I S(x(r^C) -x^expj-rtg} V QtW <^ . (3.53) 

J C+ J Z a '■ T„ 
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The integral over can be done as before since, by the parametrization, Ci T x a ) = 0- 
Use (3.34) in [1] to replace the point xf with x B to get 



v {p)t \ f -, M f |x s — x a | 2 1 ± 
K d (x a ,x B ) = J c l (r±)(w+2)/2 exp -7T rfr Xa 



= < 



( T ±)(n+2)/2 

ATg|l|x B -x a |-"l for s = 1 

J\fI(x a ,x B ; n + 2)1 for s = i 

(3.54) 

Find A/" by normalizing Kg . For s = 1, integrating over dM n and setting the result 
equal to unity gives 

B 2 -x 2 

M= ^nr- (3 - 55) 

Consequently obtain the well known result, 7 

K ° {Xa > XB) -2^R\Rn XB -xy- (3 ' 56) 



Now, for K^ the parametrization for interior points becomes 

x 2 (r,z) = [x cr (r)+ar)] 2 (3.57) 
where x cr is given by (3.6). Do the Z% integral first to get, 



f /•T ± f ~ 7r l x q'~ x n| 2 1 

Klj D \x a ,x a ,)= * a l(sr a ,)- n/2 e\ Ur a ,Vr. 

JT+ JO 



(3.58) 



Now restrict to the s — 1, n — 3 case and utilize the decomposition from Subsection 
3.2.1 in [l]. 8 In spherical coordinates, the transformation a : r a — > «r a where a = 
_R/r a takes — > 0. 9 Moreover, there is only one critical path since r a ^ excludes 
the path that passes throught the origin and intersects the boundary transversally. 
This yields 



7 Incidently, this example shows that, for the n-ball, (r^) T + = (R 2 — x 2 )/2 . 

8 For comparison purposes, this is more convenient than solving (3.58) directly (which 
entails finding a suitable normalization constant J\f after localizing the integral over T+). 

9 The transformation a must be multiplicative here and not additive as it was for the planar 
case: Additive would lead to r a > < in general, but r a t > by assumption. 
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F u (i?n Xo ,(]x a |/ J R)x a O 



r -T,\(\x a \/R)x a/ -Rn Xa \ 2 \ 



dr a 



R 



- - Rn x 



(3.59) 



The final result is 



^U D) (Xa,X a /) 



l x o| x a' 

R 



fin, 



(3.60) 



It is interesting to compare and contrast this strategy to the method of images. 

For the exterior of the n-ball, the calculation for the elementary kernel goes through as 
before. However, the boundary kernel changes sign because the normalization factor 
is found to be J\f = (x^ — R 2 )/2R since |x a | > R in this case. 



3. 3. 2 Topological n-Ball 

The procedure employed in the previous subsection can be used to determine the 
elementary kernel for more general boundaries in W 1 that are topologically equivalent 
to B n . Quite generally 10 , 



{f ) (x a ,x a ,)= / /S^W^^l^r 

J± a JO 

POO , f -^|X Q /-Xq| 2 ) 

= / l(r a 0- n/2 el ^' Ur a , 
Jo 

-/ / l(v)-" /2 el v Ur a ,VT 



(3.61) 



One recognizes the second term as a homogeneous term which ensures that the kernel 
K^ > \x a ,x a ') has the required boundary condition. 

As in the previous subsection, make a transformation(s) 11 a : x a — > x^. Then — > 
and end up with an expression that is independent of r^; 



10 I only exhibit the expressions for the case s = 1, but the general case can be handled 
analogously. 

11 Recall that one should sum over all r^ a if there are multiple critical paths. Consequently 
there may be multiple tansformations as well. 
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K { v D \x a ,x a ,) = / o °° lfa)""' 2 expj Xa| I rfr a , 

-y o !(v) 7 expj — Sdr a /. (3.62) 

Of course, the transformation(s) a can be quite complicated in general and numerical 
methods may be required, but the solution can be obtained in principle. 



3.4 The Quadrant 



As a final example, consider the upper right quadrant in M 2 , and calculate K^ 3 \x a , x a /) 
and Kf\ x a ,x a i) for the elliptic case with s — 1. This problem can be handled by 
the technique of Subsection 3.2, but it is instructive to solve it directly because it 
illustrates the case when there are multiple critical paths. 

Denote a point x G R 2 by x = (x^x 2 ). The elementary kernel is given by 



JK'{f ) (x ,x /)= / / 6{x(T a i,z),Xa>) exp{-S{x{T a ,,z))} dr a , 
Jn Uo 

6(x(T a >, z),x a ,) exp {-S(x(T a >,z))}dr a > 



(3.63) 



The first integral was already calculated in (3.11); —2 In |((x^, — x^), (x 2 , — x 2 ))|l. To 
calculate the second integral, note that there are three critical paths: a perpendicular 
line from x a to the perpendicular line from x a to the x 2 axis; and, because 

the tangent along the boundary vanishes at the origin, a straight line from x a to the 
origin. 

The transformations that take x a to each x^ are: <j\ : (x^,xl) — > (0,x 2 ); a 2 : 
(x l a1 x 2 a ) — ■> (x^,0)] and a 3 : (x^,xl) — > (0,0). The hessian of S(x cr ) is degenerate 
at the points (0,x 2 ,), (x^,,0), and (0,0). Consequently, care must be exercised when 
evaluating Fu(<r(x a ), cr _1 (x a /)). In fact the hessian vanishes at (0,0), which is there- 
fore a point of nullity 2, and one might guess that the critical path intersecting the 
origin will pick up a phase exp(i7r). 



Instead of going into details, it suffices to indicate why this expectation is reasonable 
and then simply check that it gives the answer with the correct boundary conditions. 
The integral to be evaluated is 
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Fu(0,x a / +x a ) = 




(3.64) 



Considered as an integral over the complex r a > plane, this is just iiCu(x a ,x a /) with 
x a rotated through 9 = ±7r; indicating that Fu should pick up a phase of exp(=pi7r) 
relative to Ku due to the (v) _1 factor in the integrand. 

Collecting the four contributions to (3.63) yields 



^(xa^,) = -2 {in !((*£, - x*), - **))\ - In + x'), (x 2 , - x 2 ))| 

- In |((xi, - x'), (4 + x 2 ))| + In | ((4 + xi), (x 2 , + x 2 ))|} 1 . 

(3.65) 

The first term satisfies the inhomogeneous equation, the next three terms verify the 



homogeneous equation, and K$ ^(x a , x a /) 



xi=0 



k\j \x a , x a / / 



x^=0 



0. 



Now, for the boundary kernel there are three corresponding to the three critical 
paths. For each critical time, it is best to expand about the associated critical path. 
For example, expand about 



^ r (r)=xi(l-r/r x - L J 
[4 r {r)=y? a {l-r/Ti)+xl,{r/ri) 



(3.66) 



for the critical time from x a to the x 1 axis. Each critical path will contribute to a 
segment of the boundary. Hence, the path integral separates into three terms. The 
integral has already been evaluated in (3.54) (for a different critical path). Use the 
three parametrizations for the critical paths to get 



Kf\^ o)) =M7r- i [(xi, - x\f + (xirr 1 
Kf\* a , (o,^)) =M 2 ^[{xif + {xi - xirr 1 

Kf(x a ,(0,0))=M7r- 1 [^) 2 + ^) 2 ]- 1 • (3.67) 
Finally, integrate each over its relevant range to normalize K<°>; 

^ = 2|ta,rHxVxy + V2] ; JWi= 2[ ta n-'(4/xi)+V2] ! V " = ° ' <3 - 68) 



Similar reasoning can be applied to the wedge problem where the angle is no longer 
restricted to n/2. 
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4 Conclusion 



Specialization of the general path integral constructed in [1] leads to path integral 
solutions of elliptic, parabolic, and hyperbolic PDEs with Dirichlet /Neumann bound- 
ary conditions. The path integral, along with some calculational techniques, was used 
to evaluate selected kernels of some known planar and spherical boundary problems 
in M. n . The techniques for calculating the kernels for the n-ball and quotient spaces 
are new to my knowledge, and they seem to be useful for evaluating kernels with 
more complicated geometries. It may be useful to try to solve for the kernels with- 
out recourse to the techniques developed here by applying a perturbation technique 
to . Non- vanishing V(x) can also be included and the resulting integrals can be 
solved exactly in select cases or by expansion techniques. It would be useful to do 
more substantial calculations of kernels on non-trivial manifolds. This would require 
parametrizing the paths using the Cartan development map and, probably, new cal- 
culational techniques. 



Acknowledgment I thank C. DeWitt-Morette for helpful suggestions and discus- 
sions. 



A Exponential Map Parametrization 

There are two maps that are typically used to parametrize paths on a manifold; 
the Cartan development map and the exponential map. The Cartan development 
map is presented in detail in [20]. The exponential map is particularly useful for 
parametrizing paths on linear spaces and compact Lie groups. 

The exponential map, Exp : T Xf PjM — > PjM, is a frequently used method of 
parametrization that yields the usual loop expansion when TM is a linear space. 
It maps the tangent space of the space of paths PjM at some fiducial path Xf to the 
space of L 2 ' 1 paths on M with fixed end-point at t a . 

Assume the tangent space T Xf PjM is generated from paths Xf(X, t) representing one- 
parameter variations of paths in M with fixed end-point at t a . Each variation yields 
a vector field £ = dxf(\)/d\\\ =0 along £/(t). The vector £ e T Xf PjM is a map 
C Xf : T -> T X/(t) M such that C X/ (t ) = 0. ' 

The exponential map Exp is defined point-wise on the tangent bundle in terms of the 
map exp : T X ^M — > M. For each C Xf (^) e T Xf ^)M, the corresponding family of paths 
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x(t, C; A) G M is determined by the differential equation 



dx(t; A) 
dX 



A=0 



C Xf (x(t)) x(t;0) = x f (t) 



(A.l) 



with the additional requirement x(t a ; A) = x/(t ) VA G [0, 1]. The solution in an open 
set U C M can be formally written as 



Of course, the global uniqueness and existence of the solution depends on the nature 
of TU C TM. This is a particularly useful parametrization when M is a compact Lie 
group. 

If TU is a linear space, then the parametrization reduces to 



This is the parametrization used in the examples of Section 3 with A = 1 and Xf — 
x cr since TU = R 2ra . If A is a small parameter, then one replaces x(t, £; A) with 
Xf(t) + AC X/ (t) in the path integral and expands in powers of A which leads to a loop 
expansion in the path integral. 



B Quotient Spaces 

When the manifold M can be contructed as the base space of a principal fiber bundle 
with a discrete structure group, the kernels on M can be expressed in terms of the 
kernels on the fiber bundle. The first realization of this fact came in [21] with a 
treatment of multiply connected spaces. Their result can be trivially extended to 
quotient spaces with a principal bundle structure. The construction is sketched in [22] 
and [6] : The details are included here for completeness. For a system with symmetry, 
the result allows the kernels on the base space to be expressed in terms of kernels on 
the principle bundle. In particular, using this result, it is possible to express kernels on 
multiply connected spaces, orbifolds, compact Lie groups, and homogeneous spaces 
in terms of kernels on associated covering spaces. 

Let the manifold M be a principal fiber bundle endowed with a connection, a pro- 
jection II : M -> M, and a structure group G. (If G = 7Ti(M) this yields the mul- 
tiply connected case.) Given a set of open coverings {Uj} choose a trivialization 
{Uj,^} of M and the canonical section Sj = Tp^ 1 o Id where Id : — > Uj x G and 
Tp i : n _1 (Uj) — > Uj x G. Then, under the trivialization, x = s(x)-g i— > (x, g) for x G M, 
x G M, and g G G. 




(A.2) 



x(t,C;A) =x / (t) + AC Xf (t) . 



(A.3) 
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Now let T(p )9 )M denote the tensor bundle weakly associated to M. With a choice of 
trivialization s« = (p^old on T( pg )M, define an equivariant r-form & : M — > (C m )P+<? 

#(x<7) = pQr 1 ) #(x) = pig- 1 )*- 1 ^) =■■ p{g- x )p{g) * (x) (B.i) 

for some g G G and fixed g <E G where V : M -> T (Pi?) M, # : M -> (C m ) p+9 , 
and p : G — > GL((C m ) p ^ 9 ) is a (possibly non-faithful) representation of G. Note 
that x is both a point in M and an admissible map x : (C m ) p+9 ) — > 7r _1 (M) where 7r : 
T( Pi g)M — > M. Also, p(g) \P is, by definition, the tensor field p(g) \P := x _1 i/> associated 
with the equivariant map \P. The p(g) factor has been included because there is no 
a priori relationship between 3*(x) and fl*(x). The elementary and boundary kernels 
are required to be equivariant in their arguments in the above sense. 

Consider the paths x G -PjU parametrized by the composition P = Pu ° P where 

P : Z* -> U and P u : PjU -> PjU where U C M and U C M. For a fixed 
initial point x a , a path x(-,z) can be horizontally lifted by the connection yielding 
the family of paths {5(p)(-, %>))} indexed by p where p runs over the elements g p of 
G. The (multi-) index p is discrete or continuous depending on the nature of G. The 
paths are related by £(«,)(•, %>)) = S (p) (-, z (p) ) • 5f (gp) for some g (qp) G G. 

Define the evaluation map e : Z% — > U and its lift £ : — > U. Then £ = II -1 o £. 
Paths with fixed end-point x a in the base space are lifted to paths with fixed end- 
points in the fiber over x a . It is clear that the union of the inverse images of the 
lifted paths £ _1 (x( p )(-, z^) is mapped to base space paths with fixed end-point x a by 

the evaluation map e. Hence, a path integral over PjU (which is defined in terms 
of Z%) decomposes as the sum (or integral), over the elements of a discrete (or Lie) 
group G, of integrals over Z% p where 5( p )(-, Z( p )) G . Under this decomposition, the 

elementary kernel on M becomes 




x exp {-S(x(p)(T a >, Z(p)))} dr a > VQ, p (B.2) 

where the summation/integral symbol means J2 P f° r discrete or disconnected G and 
J G dg for G the connected component of a Lie group. 

For some fiducial path with end-points X( ) a G n _1 (x a ) and X( ) o ' G n _1 (x a /), then 
X(p)(T a ',Z(p)) = x ( o)a • Z(T a >,Z(p)) • g p and x a / = x (0 ) a ' • g p ■ This, together with the 
requirement that K be equivariant in its two arguments, implies 

5(x(T a i,Z{p)),yL al ) =<5(x( )a • g v - r p (T a /,Zo),X(o)a' 

= P(^p 1 )*( X (0)o • ^p(V,^o),X( )a') (B.3) 
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where 

£ p (t ', Zq) := g~ 1 • r(v, 2q) • g p = Efa, ^ (p) ) . (B.4) 

Hence, 



K v (x a ,x a ,) = 2j)(g p 1 ) / ^(r^ -T o /)(5(x (0 ) a • r p (r a /,zo),x ( o) a ') 

*j^p «/ Op «/ C-|- 

x exp |-<S(5(p)(r a /, -%)))} dr a / £>fi p . (B.5) 
Using (B.l), the kernel on M can be expressed as 

K v (x a ,x al ) = p(^ 1 )if u (x a ,x a /) 

= P(F X ) Tfigp 1 ) I I 0{T^ - T a/ )S(x {0 )a ■ Zp(T a >, 2o),X (0 )aO 
«/p J \lp J 0-|- 

x exp |-«S(5 (?) )(T a /, ^(p)))} c?r a / £>fi p 

= : £o(^ 1 )p(F 1 )^U(p)(X(0)a,X (0 )a') 

= &(9p 1 ) K V(p)(x(0)a,X(0)a>) (B.6) 

for some g e G. In the last line, the pig' 1 ) was absorbed in a redefinition of the 
fiducial end-point (which was arbitrary anyway). 

Equation (B.6) is still not satisfactory because the term p(g p l ) is not explicitly 
known. 12 Also, there is ambiguity in the choice of the fiducial path so I could just as 
well replace g p by, say, ~g p = g v ■ g for some g G G. Proceed by taking the absolute 
value of (B.6) putting g p g p ; 



K v (x a ,x a ,) 



'^J ) (9 P 1 )-K"u(p)(x(o)o,X(o) ') 



(B.7) 



Since the modulus of the kernel should not depend on the choice of fiducial path, 
then | IT | = |itf|. A necessary and sufficient condition for this to hold is that p must 
furnish a unitary representation of G (see e.g. [21]). 

Because the fiducial path is arbitrary, it is convenient to choose X( ) = (x, e) where e 
is the identity element in G. With this choice, the elementary kernel on U is given by 

K r v (x a ,x a ,) = ^p r (g^)K v(p) (x a ,x a ,) (B.8) 



12 However, in quantum mechanical applications where & represents a wave function, (B.l) 
implies immediately that p must be a unitary (projective)representation. 
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where the p r furnish unitary representations of the structure group G and I am 
justified in writing 

JST u (x ,x /)= / / 6>(t^ - T a/ )8(x a ■ E p (T a ,,Z ),X a/ ) 

J S7p J C_|_ 

x exp {-«S(z(p)(T ',2(p)))} dT a ,VVt v . (B.9) 
where reference to the point e has been omitted since it is inconsequential. 
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